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O ! Abstract. We use uhraviolet space-based {FUSE, HST) and optical/IR ground-based (2.3m MSSSO, NTT) 

spectroscopy to determine the physical parameters of six WC4-type Wolf-Rayet stars in the Large Magellanic 
Cloud. Stellar parameters are revised significantly relative to Grafener et al. (1998) based on improved observations 
r*^ . and more sophisticated model atmosphere codes, which account for line blanketing and clumping. We find that 

CIh' stellar luminosities are revised upwards by up to 0.4 dex, with surface abundances spanning a lower range of 0.1< 

Q ! C/He < 0.35 (20-45% carbon by mass) and O/He<0.06 (<10% oxygen by mass). 

t I ■ Relative to Galactic WC5-8 stars at known distance, and analysed in a similar manner, LMC WC4 stars possess 

' systematically higher stellar luminosities, ~0.2 dex lower wind densities, yet a similar range of surface chemistries. 

, We illustrate how the classification C ill A5696 line is extremely sensitive to wind density, such that this is the 

II" ' principal difi'erence between the subtype distribution of LMC and Galactic early-type WC stars. Temperature 

. 5^ I differences do play a role, but carbon abundance does not affect WC spectral types. We illustrate the effect of 

■ varying temperature and mass-loss rate on the WC spectral type for HD 32257 (WC4, LMC) and HD 156385 

' (WC7, Galaxy) which possess similar abundances and luminosities. 

Using the latest evolutionary models, pre-supernova stellar masses in the range 11-19 Mq are anticipated for 
LMC WC4 stars, with 7-14 Mq for Galactic WC stars with known distances. These values are consistent with 
pre-cursors of bright Type-Ic supernovae such as SN 1998bw (alias GRB 980425) for which a minimum total mass 
of C and O of 14Mq has been independently derived. 

Key words, stars: Wolf-Rayet - stars: fundamental parameters - stars: evolution - stars: abundances - galaxies: 
Magellanic Clouds 



1. Introduction 

The ultimate fate of the most massive stars is likely to be 
a Type lb or Ic Supernova (SN) explosion. Type Ib's, con- 
spicuous for the absence of hydrogen in their spectra, likely 
correspond to WN-type Wolf-Rayet (WR) stars, whilst 
WC or WO stars are thought to be responsible for Type 
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Ic SN, since both hydrogen and helium are absent from 
their (early) spectra. Are the properties of WC stars im- 
mediately prior to their proposed SN explosion consistent 
with that of such SN? Models of WR stars and SN are 
now sufficiently advanced to facilitate such comparisons 
for the first time. 

Over the past couple of decades, spectroscopic tools for 
the quantitative analysis of WR stars have advanced suf- 
ficiently (e.g. Hillier & Miller 1998) to permit the reliable 
determination of abundances (Herald et al. 2001), masses 
(De Marco et al. 2000) and ionizing fluxes (Crowther et 
al. 1999). Armed with these tools, we are now in an un- 
precedented position to investigate such stars in the Local 
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Table 1. Observing log for LMC WC4 stars, including 
narrow-band photometry from Torres-Dodgen & Massey 
(1988), and catalogue numbers fr-om Breysacher (1981, Br) 
and Breysacher et al. (1999, BAT). 
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Group and beyond (e.g. Smartt et al. 2001). Already, 
high quality observations spanning UV to IR wavelengths, 
are possible for individual stars in our Galaxy or the 
Magellanic Clouds with new instruments, such as the 
Far- Ultraviolet Spectroscopic Explorer (FUSE, Moos et al. 
2000). 

WR stars in the LMC have been the focus of sev- 
eral spectroscopic investigations. WN stars have been 
studied by Conti & Massey (1989), whilst Hamann & 
Koesterke (1998), and references therein, have determined 
their quantitative properties, for which a negligible metal- 
licity effect was remarked upon relative to Galactic coun- 
terparts (the heavy metal content of the LMC is ^ OAZq, 
Dufour 1984). Torres et al. (1986), Smith et al. (1990) 
and Barzakos et al. (2001) have compared LMC WC stars, 
with Grafener et al. (1998) presenting studies of six WC4 
stars. In the latter work, WC4 stars were found to possess 
remarkably uniform (and high) carbon and oxygen abun- 
dances, plus a wide range of stellar luminosities. Results 
were inconclusive whether the mass- loss rates of LMC WC 
stars are similar to, or lower than. Galactic counterparts. 
We return to the study of these stars in this paper, since 
we possess improved and more extensive spectroscopy, in- 
cluding far-UV and near-IR datasets, plus better mod- 
eUing tools (HiUier & Miller 1998). 

We present our new observations in Sect. 2, and discuss 
the present set of model calculations in Sect. 3. Individual 
results are presented in Sect. 4 and compared with the 
previous study of Grafener et al. Quantitative compar- 
isons are made with Galactic counterparts that have been 
analysed in a similar manner (Dessart et al. 2000; HiUier 
& Miller 1999) in Sect. 5, together with evolutionary ex- 
pectations. Finally, the possibility that WC stars provide 
the precursors to Type Ic SN are discussed in Sect. 6, 
via the comparison of accurate WC masses with those of 
the CO-cores determined for recent SN Ic explosions (e.g. 
Iwamoto et al. 1998, 2000). 

2. Observations 

For our study, we utilise archival UV Hubble Space 
Telescope (HST) datasets, together with previously un- 
pubhshed far-UV FUSE, optical 2.3m Mt Stromlo & 
Siding Spring Observatory (MSSSO) spectroscopy, plus 



near-IR European Southern Observatory (ESO) New 
Technology Telescope (NTT) observations, discussed in 
turn below. An observing log for each of the six program 
stars, together with photometry from Torres-Dodgen & 
Massey (1988), is presented in Table Spectral types are 
WC4 in all cases, following either Smith et al. (1990) or 
Crowther et al. (1998). 

2.1. Far UV spectroscopy 

HD 32402, HD 37026, and HD 37680 were observed as 
part of the FUSE Principal Investigator team program 
to study hot stars (P117, P.I.: J. B. Hutchings) between 
2000 February-2001 November, as indicated in Table ^ 
The spectra were obtained as time-tag mode integrations 
of 4.5-8.1 ks duration through the 30" x 30" (LWRS) aper- 
ture. As described by Moos et al. (2000) and Sahnow et al. 
(2000), FUSE daXa consist of spectra with spectral resolu- 
tion of ^15,000 from two Lithium Fluoride (LiF) channels, 
which cover AA 990-1187 A and two Silicon Carbide (SiC) 
channels, which cover AA 905-1105 A. Spectra from each 
channel were processed by the current version of the stan- 
dard calibration pipeline (calfuse 2.0.5), which corrects 
for drifts and distortions in the readout electronics of the 
detectors, removes the effects of thermally induced grating 
motions, subtracts a background image, corrects for resid- 
ual astigmatism in the spectrograph optics, and applies 
flux and wavelength calibrations to the extracted spectra. 
Spectra from the individual channels were subsequently 
aligned, merged, and resampled to a constant wavelength 
step of 0.1 A in the manner described by Walborn et al. 
(2002b). 

2.2. Near UV spectroscopy 

All program stars were observed with the HST Faint 
Object Spectrograph (FOS) instrument between 1994 
September-1995 June (PI: D.J. HiUier). Individual expo- 
sures were obtained with the G130H, G190H and G270H 
gratings, respectively, providing complete spectral cover- 
age in the AAl 140-3301 range. These data have previously 
been discussed and published in part by Grafener et al. 
(1998). 

2.3. Optical and far-red spectroscopy 

We have used the Double Beam Spectrograph (DBS) at 
the 2.3m MSSSO telescope to observe our six targets be- 
tween 1997 Dec 24-27. Nearly complete coverage in the 
optical and far-red region was obtained in two separate 
exposures, which were in the range 240-1500 sec. The 
dichroic together with the 3001/mm (blue) and 3161/mm 
(red) gratings provided spectroscopy of 3620-6085A and 
6410-8770A, whilst the 6001/mm (blue) and 316 1/mm 
(red) gratings permitted observations at 3240-4480A and 
8640-1 lOlOA simultaneously. The detectors for both arms 
of the DBS were 1752x532 pixel SITE CCD's. A 2" slit 
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Fig. 1. Panels (a-c): Synthetic spectra for selected Civ line profiles with parameters ~ 90kK, logL/Lo=5.70, 
Uoo=3000 kms~^, C/He=0.14, O/He=0.02 by number. An unclumpcd model is shown in solid (/=!, with M = 
1.1 X IO^^Mq yr"^), moderately clumped as a dotted line (/=0.1, with M — 3.5 x 10~^Mq yr"-'^)), and highly 
clumped as a dashed-line (/ = 0.01, with M = 1.1 x 1O~^M0 yr"-"^. Electron scattering wings, and line intensities 
for the case of AA5801-12, are very sensitive to clumping, as is the shape and location of AA1548-51. Panels (d-f): 
Comparison of moderately clumped (/=:0.1) models with observations of HD 32402, showing generally good agreement. 
Ordinate units are ergcm^'^ s^^ A~^. 



provided a 2 pixel spectral resolution of ^SA, whilst a 6" 
slit (for absolute spectrophotometry) provided a spectral 
resolution of ^12 A. A standard data reduction was car- 
ried out, including absolute flux calibration, using wide slit 
spectrophotometry of HD 60753 (B3IV) and ^ Col (09 V), 
plus atmospheric correction using B stars HR2221, 4074 
and 4942. 

2.4. Near-IR spectroscopy 

Long slit, near-IR spectroscopy of HD 32402 was acquired 
on 1999 Sept 1-2 with the NTT, using the Son of Isaac 
(SofI) instrument, a 1024x1024 pixel NICMOS detector, 
and low resolution IJ (GRB) grating, with spectral cover- 
age 0.94-1.65 /im and a dispersion of 7.0 A/pix. The 0.6" 
slit provided a 2-pixel spectral resolution of 14A. The total 
integration time was 480 sec (1 Sept) and 960 sec (2 Sept). 
The removal of atmospheric features was achieved by ob- 
serving HD 38150 (F8V) immediately after HD 32402, 
at a close airmass (within 0.03). Similar observations of 
Hip 22663 (B2 V) permitted a relative flux correction, us- 
ing a Tcff =22kK Kurucz model atmosphere normalized to 
V=7.66 mag. 

A standard extraction and wavelength calibration 
was carried out with iraf, while FIGARO (Shortridge, 



Meyerdierks & Currie, 1999) and dipso (Howarth et al. 
1998) were used for the atmospheric and flux calibration, 
first artificially removing stellar hydrogen features from 
the HD 38150 spectrum. Our relatively fluxed dataset was 
adjusted to match MSSSO spectrophotometry. 

3. Model Analysis 

3.1. Method 

For this study, we employed the code of Hillier & Miller 
(1998), CMFGEN, which iteratively solves the transfer 
equation in the co-moving frame subject to statistical and 
radiative equilibria in an expanding, spherically symmet- 
ric and steady-state atmosphere. These models account for 
clumping, via a volume filling factor, /, and line blanket- 
ing, since these effects affect significantly the derived prop- 
erties of WC stars (e.g. ffilHer & Miller 1999). Through the 
use of 'super- levels', extremely complex atoms can been 
included. For the present application, a total of 2526 lev- 
els (combined into 788 super- levels), 40 depth points and 
26 239 spectral lines of He, C, O, Ne, Si, S, P, Ar and 
Fe are considered as indicated in Table ||. The reader is 
referred to Dessart et al. (2000) for details of the oscilla- 
tor strengths, collision and photoionization cross-sections. 
The OPACITY (Seaton 1987, 1995) and IRON projects 
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Table 2. Summary of the WC model atom calculations. N/ is the number of full levels, the number of super levels 
and Nt the corresponding number of transitions. A total of 2526 levels (788 super levels) and 26 239 transitions are 
considered. 
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(Hummer et al. 1993) formed the basis for most atomic 
data utilised in this work. Test calculations, addition- 
ally involving extensive model atoms for Naii-vil, Mgiii- 
VII, CI IV- VII, Caiv-ix, Criv-vi, Mniv-vii and Niiv-ix 
were also considered, with minor effects on the emergent 
spectrum, except for an increase in the strength of Civ 
AA1548-51, 5801-12. 

Following theoretical discussions by Schmutz (1997), 
and for consistency with previous studies, we adopt a form 
for the velocity law (Eqn 8 from Hillier & Miller 1999) 
such that two exponents are considered (/3i = l, /32=20), 
together with an intermediate (uext) and terminal Uoo ve- 
locity, with Uoxt ~ Voo - 300 kms""'^. Consequently, acceler- 
ation is modest at small radii, and continues to relatively 
large distances, i.e. 0.9uoo is reached at ^20i?* versus lOi?* 
for a standard P—1 law. In the future, we intend to re- 
lax this assumption via hydrodynamical driving, at least 
for the outer envelope. Terminal velocities are measured 
principally from the C iv AA1548-51 UV P Cygni profiles, 
subject to minor adjustment at the line fitting stage, and 
agree relatively well (within 15%) with previous determi- 
nations from Grafener et al. (1998) who measured wind 
velocities using the unsaturated C ill A2297 line. 

The formal solution of the radiative transfer equation, 
yielding the final emergent spectrum, is computed sepa- 
rately. Generally, Doppler profiles are adopted through- 
out, except for C ill A977, C iv AA1548-51 for which Voigt 
profiles are absolutely essential. Except where noted, these 
final calculations assume uturb^SO kms~^. 

A series of models were calculated for each star in 
which stellar parameters, T*, log(L/LQ), M/^/f, C/He 
and 0/He, were adjusted until the observed line strengths 
and spectroscopic fiuxes were reproduced. Our spectro- 
scopic analysis derives M/y/J, rather than M and / in- 
dividually, since line blending is so severe in WC4 stars. 
This is illustrated in Fig. |^(a-c), where we compare se- 
lected synthetic line profiles using identical models of 
HD 32402, except that volume fiUing factors of /=1.0, 0.1 
and 0.01 are adopted (with corresponding changes in M). 
C IV A7724 is typical of the vast majority of lines, in that 
solely changes in electron scattering wings are predicted. 



Emission intensities of some lines, such as AA5801-12, are 
very sensitive to clumping, as are the shapes of the car- 
bon resonance lines CiV AA1548-51 and Cm A977. Recall 
that a Voigt profile is necessary for these lines. It is appar- 
ent from Fig. 0(d-f) that the /=0.1 case provides the best 
overall match to observed electron scattering wings, and 
to the central wavelength of AA1548-51, which unusually 
distinguishes itself as a clumping probe in WC stars in 
this way. 

Hamann & Koesterke (1998) considered clumped mod- 
els of HD 37026, with /=1.0, 0.25 and 0.0625, and came 
to the conclusion that the latter provided the best fit to 
observations. In general, we can firmly rule out /=1.0 
for all cases, but our (admittedly simplistic) treatment 
of wind clumping prevents a precise determination - see 
Dessart & Owocki (2002ab) for recent progress in this 
field. Ultimately, we cannot pretend to derive the dumpi- 
ness of the wind, and so we have selected identical values 
of /(=0.1) for the other program stars, as a reasonable 
compromise to the observed spectra in each case, and for 
consistency with other (Galactic) stars analysed in our 
previous studies. 

3.2. Stellar Parameters and Abundances 

The wind ionization balance is ideally deduced using iso- 
lated spectral lines from adjacent ionization stages of 
helium, carbon or oxygen. In order to ensure consis- 
tency with studies that are restricted to optical or far-red 
wavelengths we selected Hei 10830/Heii 5412 and Cm 
6740-h8500-h9710/Civ 7700. No suitable oxygen diagnos- 
tics were available in the optical or far-red, whilst use of 
He I 10830 was difficult in practise for our program stars 
(given the poor CCD response in the near-IR), with the 
exception of NTT-Sofl data for HD 32402. Consequently, 
carbon provided our principal diagnostic lines, with the 
carbon-to-helium abundance determined from C iv 5471 
and Hen 5412. We consider the accuracy of our derived 
luminosities to be 15%, with effective temperatures some- 
what worse constrained. Mass-loss rates also have a 15% 
accuracy if the adopted filling factor if correct. 
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Fig. 2. Comparison between the de 
troscopy of HD 32125 (solid) in the He 
5471 region, and theoretical predictions 
dition to those features identified {W\ 
numerous weaker carbon, oxygen, neon 
tures. From top-to-bottom: C/He=0.07 
by mass), 0.13 (25%, adopted) and 0.21 
units are ergcm"^ s""'^ A~^. 
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In Fig. |2| we present de-reddened observations of 
HD 32125, plus three identical models except that C/He 
ratios are 0.07, 0.13 (adopted) and 0.21 by number, corre- 
sponding to 15%, 25% and 35% carbon by mass, with the 
oxygen mass fraction maintained at 10%. The relatively 
small changes in line strength for differing C/He ratios 
emphasise the need for high S/N and absolute flux cal- 
ibration. Difficulties are only encountered for stars with 
very large wind velocities in which Civ 5471 and Hell 
5412 are severely blended. Overall, we consider an (inter- 
nal) accuracy of 10% in our derived carbon mass fraction, 
i.e. 25±3% for the case illustrated above. External errors 
are difficult to quantify, but from test calculations carried 
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Fig. 3. Comparison between identical synthetic spec- 
tra for HD 32402, except that iron abundances are cho- 
sen to be 0.2^0 (dashed, blue), OAZq (solid, black) 
and 0.8^0 (dotted, red). This illustrates how difficult 
it is to accurately determine iron abundances based on 
UV spectroscopy in WC stars. Differences in iron con- 
tent within this range has little effect on the emergent 
carbon-oxygen-helium (diagnostic) lines. Ordinate units 
are ergcm"^ s~^ A~^. 



out using the line blanketed code of Grafener et al. (2002) 
the degree of consistency is excellent. 

In contrast with carbon, the determination of oxygen 
abundances are much more difficult. This is due to the 
absence of recombination lines amenable to simple anal- 
ysis, a much more complex ionization balance structure 
than He and C, and more problematic line blending - 
for example the well known oxygen feature between C iv 
A5471 and C ill A5696 is almost equally split between O ill 
A5592 and O V AA5572-5607. Nevertheless, there are nu- 
merous UV and optical lines of oxygen which are sensi- 
tive to abundance variations, which are discussed below 
for HD 37026. For oxygen we admit an (again internal) 
accuracy of 50% in our derived oxygen mass fractions. 
Absolutely, a factor of two is probably a more objective 
value. Again, test calculations by Grafener (priv. comm.) 
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indicates good agreement, except for Oiii lines for which 
dielectronic contributions to lines differ somewhat. 

For elemental species other than He, C and O, we gen- 
erally adopt OAZq abundances (Si, S, P, Ar and Fe). We 
have considered models in which the iron abundance is 
doubled or halved, without particularly affecting the emer- 
gent spectrum. This is illustrated in Fig. ^ for UV syn- 
thetic spectra of HD 32402. For the case of neon, there is 
observational evidence that Ne is enriched in Galactic WC 
stars such that its mass fraction is ^--^1.5% (Ne/He=0.004 
by number, Dessart et al. 2000). Consequently, we have 
sought UV/optical line diagnostics as evidence for Ne en- 
richment in LMC WC4 stars (see below). 

3.3. Reddening and Distance 

Simultaneously with the determination of the ionization 
balance, wind density and chemistry of each star, the 
luminosity is adjusted to ensure consistency with de- 
reddened spectroscopy. In all Galactic foreground 
extinction of i?B-v=0.07 is adopted following Schlegel et 
al. (1998), using the parameterization of Seaton (1979), 
with R=3.1=Ay / Eb-y . The extinction law followed by 
Seaton differs substantially from Cardelli et al. (1989) for 
A <1500A. The distance modulus to the LMC was as- 
sumed to be 18.50, as derived by Cepheids and SN 1987A 
(Lee et al. 1993; Panagia et al. 1991), corresponding to 
50kpc. 

LMC contributions are varied in -Eb-v a-nd R following 
the parameterization of Howarth (1983). In the far-UV, in- 
terstellar absorption in the H i atomic Lyman series is ac- 
counted for following the technique of Herald et al. (2001). 
Correction for molecular hydrogen in the FUSE spectral 
region has been made by using H2 column densities deter- 
mined by Tumlinson et al. (2002): log (H2/cm"2)=15.45 
and 18.94 for HD 37026 and HD 37680, respectively The 
H2 column density towards HD 32402 is substantial, with 
log (H2/cm-2) 0.19. HD 37026, therefore, has by far the 
'cleanest' far-UV stellar spectrum of the three WC4 stars 
thus far observed with FUSE. 

4. Results for individual stars 

We have decided against presenting spectral fits for each 
star in turn, since the model successes and failures are gen- 
erally common to all our program stars. We solely present 
spectral fits to the entire far-UV, UV, optical and near-IR 
observations of one star, HD 37026 in Fig. ^. 

4.1. Study of HD 37026 (Brey 43 = BAT99-52) 

Overall, the flux distribution of HD 37026 is very well re- 
produced by the synthetic model - particularly the iron 
forest in the AA1250-1500 region and plethora of near- 
UV and optical carbon-helium features with C/He=0.32 
by number (44% by mass). Hei A10830 may be rather 
underestimated in strength, although the low S/N of the 
observations, due to the poor CCD response beyond 1/im, 



hinders detailed comparisons. To illustrate the reliability 
of ncar-IR synthetic spectra for WC4 stars based on stud- 
ies from optical diagnostics we present high S/N NTT-SofI 
observations of HD 32402 in Fig. || with model predictions. 
Overall the match is excellent, including Hei A10830. 




I ■ ■ • ■ ■ ^ 

1.5 1.8 2.1 

wavelength (/xm) 

Fig. 5. Comparison between de-reddened near-IR NTT- 
Sofl observations of HD 32402 (solid, black) and synthetic 
spectra (dotted, red) determined from optical diagnostics. 
Ordinate units are ergcm"^ s^^ A~^. 

Problematic carbon features for HD 37026 presented in 
Fig. ||are as follows. CiV AA1548-51 emission is predicted 
to be somewhat too weak. Previous difficulties with repro- 
ducing the Civ AA5801-12 feature (e.g. Hillier & Miller 
1999) appear to have been mostly resolved by the inclu- 
sion of additional metal opacity (specifically Ne and Ar) . 
Unfortunately, C ill lines provide conflicting information. 
For example, A977, A1909, A2297, AA4647-50 and A8500 
are reproduced to within 25% in HD 37026. However, 
A6470 is up to 50% too weak, whilst A9710 is similarly too 
strong. The blend including C ill A3747 is systematically 
too weak, despite efforts to improve atomic data. Amongst 
other WC4 stars, consistency between A6470 and A9710 is 
generally excellent. 

As introduced above, there are numerous UV and 
optical lines of oxygen which are sensitive to abun- 
dances variations although results are somewhat dis- 
crepant. Fig. ^ presents identical models for HD 37026 ex- 
cept that the oxygen mass fraction is 5% (O/He=0.025), 
10% (O/He=0.05) or 20% (O/He=0.1). Diagnostics which 
correlate well with abundance are Ov A1815 (indicating 
5%), Oiii A2983 (10%), Oiv A3072 (~8%), Oiv AA3560-3 
(~8%) and Olii-hv A5590 (>20%). It is unfortunate that 
for HD 37026, the only strong oxygen feature available in 
the visual implies a substantially different oxygen abun- 
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dance from the numerous other diagnostics. This is not 
always the case - A5590 is representative of near-UV di- 
agnostics for HD 32125, HD 32257 and HD 37680. 

Great care should be taken when comparing results 
from solely optical observations, as is common for extra- 
galactic WC stars (e.g. Smartt et al. 2001), with those 
additionally including UV and far-UV observations. In 
the case of HD 37026, we adopt a mean oxygen abun- 
dance from the above hst of diagnostics, i.e. 9% by mass. 
The most significant model discrepancy is the predicted 
strength of the O VI lines - notably AA1032-38 (for those 
stars with far-UV datasets), AA38 11-34 and A5290. We 
have considered the inclusion of X-rays to maintain the 
population of high oxygen ions in the wind, but without 
much success, especially for the optical O VI lines. We are 
only able to maintain such high ions if the wind density 
were reduced, as for Sand 2 (WO, Crowther et al. 2000), 
but at the expense of the goodness of fit for other diag- 
nostics. 

Whilst the majority of spectral features are from 
He/C/O/Fe, there are signatures of several other elements 
in the spectra of WC4 stars. These are most readily ap- 
parent from HST and FUSE observations. Prominent P 
Cygni profiles of Siiv AA1393-1402, Pv AA1118-28 and 
S IV AA1062-73 features are present, and successfully re- 
produced with abundances scaled to OAZq. 

Unfortunately, spectral features of neon (generally 
Neiii-v) in WC4 stars are rather weak, even for large 
adopted abundances. This is illustrated in Fig. ^ where 
UV and optical observations of HD 37026 are compared 
with models accounting for a varying amount of neon 
- 2% (dotted), 0.4% (dashed), and 0% (dot-dash). The 
strongest accessible (unblended) UV features lie in the 
vicinity of Cm A2297 - Neiv AA2043-73, AA2204-21, 
A2175, AA2351-63, Nev A2260. These generally favour a 
neon abundance of <0.4%, which still represents a poten- 
tial overabundance by a factor of <5. Optical diagnostics 
are even harder to locate. Neiv A5045, A5245 are amongst 
the strongest features, but cannot readily be used to de- 
termine neon abundances without a superior match to the 
observed spectrum from other species. We are therefore 
unable to verify evolutionary predictions for the Ne abun- 
dance in WC stars beyond those discussed by Dessart et 
al. (2000) from mid-IR fine structure lines. 

The temperature, density, radius and ionization bal- 
ance of He, C, O, Ne and Fe for our final HD 37026 model 
versus wind velocity are presented in Fig. |^. This illus- 
trates the complex ionization structure for heavy elements 
with respect to helium and carbon. 

4.2. Summary of results and comparison with Grafener 
et al. 

Deduced stellar properties for each star are presented in 
Table ^, including ionizing fluxes below the Hi {Qo), Hei 
(Qi), and Hen {Q2) edges. We also present final synthetic 




2000 2100 2200 2300 2400 





Neiv 


Neiv 
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wavelength (A) 

Fig. 7. Comparison between de-reddened observations of 
HD 37026 (solid, black) and synthetic spectra for a variety 
of neon abundances - 2% by mass (dotted, blue), 0.4% 
(dashed, red), and 0% (dot-dashed, green). Ordinate units 
are ergcm"^ s~^ A~^. 

spectra with each observation for the diagnostic He ii 5412 
and C IV 5471 profiles in Fig. |. 

The present sample of stars has previously been stud- 
ied by Grafener et al. (1998), whilst Hamann & Koesterke 
(1998) have considered clumped models for HD 37026. 
We present a comparison of the derived parameters from 
these various studies in Table |[ Overall, we obtain sim- 
ilar stellar temperatures to Grafener et al., but substan- 
tially higher stellar luminosities (0.1-0.4 dex), lower wind 
densities (0.2±0.1 dex, after correcting for the clumped 
wind) and higher wind velocities (up to ^15%). We concur 
with Grafener et al., in that LMC WC4 stars form a rela- 
tively homogeneous sample, with high luminosities, tem- 
peratures and comparable wind densities. However, our 
improved observational and theoretical approach has led 
to substantially revised stellar parameters, notably relat- 
ing to chemical abundances. Our improved optical spec- 
troscopy permits a much finer determination of C/He ra- 
tios, e.g. for HD 37680 we obtain C/He=0.10 by number 
versus C/He=0.27 according to Grafener et al. - see Fig. ^ 

Differences between the two methods are as follows: 

1. Obscrvationally, Grafener et al. utilised identical 
HST/FOS UV spectroscopy, plus 10-15A resolution 
optical datasets from Torres & Massey (1987) cover- 
ing AA3375-7375, with the exception of HD 32402 for 
which they had higher quality ESO spectra. In all cases 
their observations omitted the useful far-red MSSSO 
and far-UV FUSE diagnostics. 

2. Theoretically, similar model assumptions were made 
by Grafener et al., except that relatively simple model 
atoms of helium, carbon and oxygen were used (e.g. 
O III was omitted) . We deduce substantially lower 
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Table 3. Stellar parameters for the program LMC WC4 stars, plus Sand 2 (WO) from Crowther et al. (2000) analysed 
in the same manner. Note that all models assume a filling factor, /, of 0.1, whilst a two-component velocity law was 
used in all cases with /3i — 1, (32 = 20, Woxt = (^oo — 300) km s~^. 



HD 


T. 


7?* 


log Lph 


T2/3 


R.2/3 


logM 




log Qo 


logQi 


logQ2 


C/He 


O/He 


/3ho 


/3c 


Po 




kK 


Rq 


Lq 


kK 


Rq 




kms~^ 


s-i 


s-i 


s-1 






% 


% 


% 


32125 


90 


2.2 


5.44 


74 


3.2 


-4.8 


2500 


49.32 


48.91 


39.03 


0.13 


0.04 


65 


25 


10 


32257 


85 


2.3 


5.42 


71 


3.4 


-4.9 


2300 


49.29 


48.81 


36.25 


0.35 


0.04 


45 


47 


8 


32402 


90 


2.9 


5.70 


67 


5.2 


-4.5 


3000 


49.55 


49.10 


36.11 


0.14 


0.02 


66 


28 


5 


37026 


90 


2.8 


5.65 


68 


4.8 


-4.5 


2900 


49.51 


48.98 


38.80 


0.32 


0.05 


46 


44 


9 


37680 


85 


3.2 


5.68 


62 


6.0 


-4.4 


3200 


49.53 


49.04 


38.36 


0.10 


0.01 


74 


22 


3 


269888 


85 


2.4 


5.44 


72 


3.4 


-4.8 


2600 


49.31 


48.86 


36.27 


0.32 


0.06 


45 


43 


11 


Sand 2 


150 


0.65 


5.28 


110 


1.2 


-4.9 


4100 


49.10 


48.85 


40.31 


0.70 


0.15 


27 


56 


16 



O/He ratios (<0.1) than Grafener et al. (0.05-0.25), 
plus a lower range of C/He ratios due to our improved 
optical spectroscopy (0.1-0.35 versus 0.3-0.55) and our 
strict adherence to He 11 5412/CiV 5471 as abundance 
diagnostics, plus the allowance for clumping. It is pos- 
sible that erroneous C/He abundances may be derived 
if the He 11 5412 electron scattering wing is unrealisti- 
cally strong. 

3. Line blanketing was neglected by Grafener et al., as 
was the contribution to the opacity by iron and other 
heavy elements. In contrast, we account for Ne, Si, P, 
S, Ar and Fe, and also consider the (minor) influence of 
other heavy elements. The neglect of blanketing, and 
use of restricted model atoms yields erroneous rectified 
synthetic spectra both in the UV (due to the iron for- 
est) and the optical (due to many broad overlapping 
features), see Hillier & Miller (1998; 1999). 

4. Iron line blanketing is principally responsible for the 
increased stellar luminosities obtained here relative 
to Grafener et al. We obtain bolometric corrections 
in the range —4.8 < B.C. < —4.5 mag, up to 0.5 
mag greater than WC6-8 stars, but 1 mag lower than 
Sand 2 (Crowther et al. 2000). We also derive some- 
what higher interstellar reddenings to the program 
stars than Grafener et al. via our spectral synthesis, 
typically by Ay ~0.2 mag as presented in Table ||, 
owing to the superior technique followed. 

Recently, Grafener et al. (2002) presented a line blan- 
keted study of HD 165763 (WC5) where many of their pre- 
vious deficiencies are rectified, i.e. enhanced model atoms 
for carbon and oxygen, line blanketing from iron group 
elements, clumping, and a more objective line profile fit- 
ting technique. Indeed, test calculations carried out by 
Grafener (priv. comm.) for HD 37026 show excellent con- 
sistency with the present results, giving confidence in re- 
sults obtained with the two independent codes. 

Whilst Grafener et al. included the opacity from other 
iron group elements, Ne . . . Ar were omitted. We find that 
additional opacity sources (namely Na, Mg, CI, Ca, Cr, 
Mn, Ni) play negligible spectroscopic role, other than a 
strengthening of Civ AA1548-51, AA5801-12 emission, but 



Table 4. Quantitative comparison between fundamental 
stellar parameters newly derived here for LMC and (se- 
lected) Galactic WC stars (referred to as C02) - using 
line blanketed, clumped models of Hillier & Miller 1998) 
- and by Grafener et al. (1998, G98) and Koesterke & 
Hamann (1995, KH95) using non-blanketed homogeneous 
models, and Hamann & Koesterke (1998 HK98) using non- 
blanketed clumped models for HD 37026. Catalogue num- 
ber for LMC (BAT: Breysacher et al. 1999) and Galactic 
(WR: van der Hucht 2001) stars are indicated in paren- 
thesis. 



HD 


log Lph 

kK Lq 


log(^) 


kms~^ 


C/He 


O/He 


D 

kpc 


Ref 






LMC 










32125 


90 5.44 


-4.3 


2500 


0.18 


0.04 


50 


C02 


(9) 


95 5.29 


-4.2 


2300 


0.32 


0.12 


50 


G98 


32257 


85 5.42 


-4.4 


2300 


0.35 


0.05 


50 


C02 


(8) 


104 5.13 


-4.1 


2300 


0.43 


0.24 


50 


G98 


32402 


85 5.70 


-4.0 


3000 


0.14 


0.02 


50 


C02 


(11) 


92 5.62 


-3.8 


2800 


0.55 


0.16 


50 


G98 


37026 


90 5.65 


-4.0 


2900 


0.32 


0.05 


50 


C02 


(52) 


103 5.26 


-3.9 


2600 


0.43 


0.24 


50 


G98 




141 5.6 


-3.7 


2800 


0.43 


0.24 


50 


HK98 


37680 


85 5.68 


-3.9 


3200 


0.10 


0.01 


50 


C02 


(61) 


94 5.55 


-3.7 


2800 


0.27 


0.05 


50 


G98 


269888 


85 5.44 


-4.3 


2600 


0.32 


0.06 


50 


C02 


(90) 


104 5.13 


-4.1 


2300 


0.43 


0.24 


50 


G98 






Galactic 










76536 


80 5.38 


-4.2 


2055 


0.25 


0.03 


2.0 


C02 


(14) 


86 4.9 


-4.2 


1800 


0.14 




2.0 


KH95 


213049 


80 5.02 


-4.5 


2280 


0.30 


0.04 


2.7 


C02 


(154) 


71 4.9 


-4.2 


2050 


0.22 




3.5 


KH95 



do contribute to the line driving force[J. From test calcula- 
tions, the inclusion of these additional lines does help, but 
the effect is still insufficient to supply the necessary driv- 

^ These additional elements are omitted from our present 
analysis for consistency with our recent studies of Galactic WC 
stars (e.g. Dessart et al. 2000). 



Crowther et al.: Study of LMC WC4 stars 



9 



Tg (10* K) 




1000 2000 3000 1000 2000 3000 




1000 2000 3000 1000 2000 3000 



6+ n5+ ^^4+ p.3 + 
























/c+\ 









0"+ 




1000 2000 3000 1000 2000 3000 




1000 2000 3000 1000 2000 3000 

V (kms"^ ) V {kms"^ } 



Fig. 8. Atmospheric structure of our final HD 37026 
model, indicating the variation of temperature (lO^K), 
density (cm^'') and radius (i?*) with velocity (kms~^), 
plus the ionization balance (in N/Ntot) for helium, car- 
bon, oxygen, neon and iron. A standard /3=1 velocity law 
is also included on the radius plot (dotted line). 

Table 5. Interstellar reddenings, absolute magnitudes, 

and H I column densities (from fits to Lya) to the pro- 
gram LMC WC4 stars, adopting a mean Galactic fore- 
ground of Eb-v = 0.07 from Schlegel et al. (1998). 
Narrow-band photometry is taken from Torrcs-Dodgen & 
Massey (1988), whilst we include comparisons with pre- 
vious reddening determinations by Grafener et al. (1998) 
and Morris et al. (1991). 
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et al. ( 
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32125 


15.02 


0.07+0.09 


0.03+0.08 


0.01 


21.5 


-4.02 


32257 


14.89 


0.07+0.08 


0.03+0.06 


0.08 


21.0 


-4.12 


32402 


13.89 


0.07+0.04 


0.03+0.05 




21.0 


-5.00 


37026 
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0.07+0.04 


0.03+0.01 


0.08 


21.0 


-4.85 


37680 


14.03 


0.07+0.06 


0.03+0.05 


0.05 


21.3 


-4.92 


269888 


15.41 


0.07+0.27 


0.03+0.28 




22.0 


-4.19 
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Fig. 9. Comparison between the de-reddened spec- 
troscopy of the targets (solid) in the Hcii 5411, Civ 
5471 and Oiii-V 5590 region, and final theoretical pre- 
dictions (dotted), demonstrating the excellent agreement 
reached for C/He determinations. Ordinate units are 
erg cm^^ s^^ A"-'^. 

Table 6. Observing log for Galactic WC stars, including 
narrow-band photometry from Torres-Dodgen & Massey 
(1988), catalogue numbers and cluster /association dis- 
tances from van der Hucht (2001). Optical spectropho- 
tometry for HD 213049 is taken from Torres & Massey 
(1987). 



HD 


WR Sp Type 


lUE 


MSSSO/ 


V d 




M„ 








KPNO 


mag kpc 


mag 


mag 


76536 


14 WC7 


Sep 80 


Dec 97 


9.42 2.0 


0.61 


-4.6 


213049 


154 WC6 


May 82 


Oct 80 


11.54 2.7 


0.72 


-3.6 



ing force, especially in initiating the outflow at the base 
of the wind, in common with Grafener et al. (2002). This 
will be discussed in a forthcoming publication. 

5. Comparison with properties of Galactic WC 

stars 

We have presented a detailed study of six LMC WC stars. 
How do their properties compare with Galactic counter- 
parts? What differences/similarities might be expected 
from evolutionary predictions? 

It is well known that WC4 stars are relatively rare in 
the Milky Way, with only 5 known out of 87 WC stars (van 
der Hucht 2001), whilst they dominate the LMC carbon 
sequence with 15 known out of 23 WC stars (Breysacher 
et al. 1999). Indeed, most LMC WC stars with WC5 6 
spectral types are WC4+0 binaries with the later appar- 
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ent subtypes resulting from C ill A5696 emission due to 
wind collision (Bartzakos et al. 2001). 

5.1. New analyses of Galactic WC stars 

It is unfortunate that only a few Galactic WC stars lie at 
known distances. There are no single cluster/association 
WC4 stars, so we are unable to directly compare our LMC 
results with Galactic counterparts. Instead, comparisons 
must be made with WC5-8 stars, which have been anal- 
ysed in a similar manner, namely 7 Vel (WRll, WC8-I-0, 
De Marco et al. 2000); HD 92908 (WR23, WC6, Smartt 
et al. 2001), HD 165763 (WRlll, WC5, HiUier & Miller 
1999), and HD 192103 (WR135, WC8, Dessart et al. 2000). 
To this short list we have added HD 156385 (WR90, WC7) 
for which Dessart et al. (2000) obtained a distance based 
on a calibration of (binary) WC7 cluster/association mem- 
bers. 

In order to produce a more substantial set of data 
for Galactic WC stars, we have newly analysed two stars 
whose distances are known from OB association member- 
ship - HD 76536 (WR 14, WC7) and HD 213049 (WR 154, 
WC6). Spectrosopic observations, cluster distances and 
derived reddenings are listed in Table ||. Ultraviolet 
datasets generally represent low dispersion International 
Ultraviolet Explorer (lUE) spectroscopy, with the excep- 
tion of three (merged) high dispersion short wavelength 
(SWP) observations of HD 76536. Optical datasets are 
drawn from Mt Stromlo 2.3m DBS observations, ob- 
tained during our LMC Dec 1997 observing campaign for 
HD 76536, plus archival Torres & Massey (1987) datasets 
for HD 213049 obtained with the 0.9m Kitt-Peak (KPNO) 
intensified Reticon scanner (IRS). We obtain somewhat 
higher reddening determinations than the recent litera- 
ture (Koesterke & Hamann 1995; Morris et al. 1993) for 
HD 76536 and HD 213049, namely £'b-v=0.42-0.49 and 
0.63-0.68, respectively. 

Identical analysis methods were employed to those pre- 
sented above. For conciseness, we omit presenting fits to 
spectral hnes for HD 76536 and HD 213049, ahhough the 
quality is comparable with those of HD 156385 (WC7, 
Dessart et al. 2000) and HD 92809 (WC6, Smartt et al. 
2001). Note that earlier problems with reproducing the 
strength of C ill A4186 have been resolved and were caused 
by a 'poor' super-level assignment. 

We present derived stellar parameters in Table ^, to- 
gether with those obtained using non-blanketed helium- 
carbon models by Koesterke & Hamann (1995). The use 
of blanketed model atmospheres result in higher stellar lu- 
minosities as has been discussed above for the LMC WC4 
stars, whilst mass-loss rates are again reduced even af- 
ter correcting for clumping. C/He ratios are increased in 
this case, although Koesterke & Hamann (1995) neglected 
oxygen in their calculations. 
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Fig. 10. Comparison between the mass-loss rates and lu- 
minosities of LMC WC4 (filled squares), WO stars (filled 
triangle, Crowther et al. 2000) and Galactic WC5-8 stars 
(open squares) at known distances. The lines correspond 
to the calibrations of Nugis & Lamers (2000) for Galactic 
WR stars (dotted), assuming C/He=0.2, C/0=4 by num- 
ber, plus a linear fit to the LMC WC and WO results 
(sohd). 

5.2. A metallicity effect? 

In Fig. |l^ we compare the luminosities and mass-loss rates 
of LMC WC stars with Galactic counterparts located at 
known distances, including HD 76536 and HD 213049 for 
the first time. Observationally, two main effects are ap- 
parent: 

— It is striking that the mean luminosity of WC stars 
in the LMC, log(L/_L0)=5.6, is a factor of two times 
higher than those of their Galactic counterparts, 
log(i/L0)=5.3, which have been studied in an analo- 
gous manner. The conventional interpretation of this 
difference is that LMC stars are descended from, on 
average, higher initial mass stars, as might be ex- 
pected from evolutionary predictions that account for 
the lower mass-loss rates during the main-sequence and 
post main-sequence phases at lower metallicities. A 
more speculative alternative is that their initial masses 
span a similar range, but that they are currently faster 
rotators owing to lower mass-loss rates during earlier 
evolutionary phases. Polarization evidence points to 
the majority of Galactic WC stars being fairly symmet- 
ric, and hence are slow rotators by inference (Harries 
et al. 1998). 

— LMC WC stars possess lower mass-loss rates than 
Galactic stars of similar luminosity, unless clumping 
factors differ between the two samples. The semi- 
empirical calibrations of Nugis & Lamers (2000) for 
Galactic WR stars match the observed WC properties 
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rather well, whilst a linear fit to the LMC WC and 
WO data reveals 



log (A/) = 1.381og(i/L0) - 12.35 



(1) 



with a similar slope to the Nugis & Lamers (2000) 
generic WR calibration (their Eqn. 22), albeit offset 
by -0.2 dex. 

The metallicity dependence of the mass-loss rate for 
O stars (M cx ^0-5~o.7^ jg .^gjj established observation- 
ally (Puis et al. 1996) and theoretically (Kudritzki & Puis 
2000; Vink et al. 2001). The situation for WN stars is 
less clear (Crowther 2000; Hamann & Koesterke 2000). 
However, since iron group elements provide most of the 
driving for the winds of WN stars, they should exhibit 
a similar dependence on metallicity as the 0-type stars. 
A metallicity dependence of — Z'^-^ would predict that 
LMC WR stars have winds that are 0.2 dex lower than 
equivalent Galactic stars, as appears to be the case. 

One might question whether a metallicity dependence 
of WC stars is expected, since their winds are so domi- 
nated by carbon and oxygen. From test calculations for 
WC4 stars, wc find that the contribution of lines from 
dominant ions (~250 eV) at the base of the wind (< 100 
kms"""^) is crucial to initiating the outflow. Carbon does 
not possess lines from such ions since the ionization po- 
tential of 0^+ (64 eV) and (392 eV) do not lie in the 
necessary range - the same is true of oxygen (recall Fig. ^ . 
Instead, ions from Ne, Ar, and Fe-group elements appear 
to be critical in initiating the outflow from early-type WC 
stars, which would qualitatively explain their sensitivity 
to heavy metal content. 

We shall see in the next section that a metallicity de- 
pendence of mass-loss rates contributes to an explanation 
for the subtype distribution amongst WC stars. In ad- 
dition, since weaker WR winds are more transparent to 
harder ionizing photons (Schmutz et al. 1992; Crowther 

1999) one would expect a greater contribution of WR stars 
to the hard ionization of young starburst regions in more 
metal deficient galaxies. The inclusion of line blanketed 
metallicity dependent mass-loss rates for WR stars in evo- 
lutionary synthesis calculations is discussed by Smith et 
al. (2002). 

The Galactic and LMC WC populations are further 
compared in Fig. |ll|, where we present their (C4-0)/IIe 
values versus luminosity, together with evolutionary pre- 
dictions for non-rotating stars (solid, Meynet et al. 1994) 
for the two metallicities, plus results for an initially rapidly 
rotating 4OM0 model at l.OZ© (dotted, Meynet & Maeder 

2000) . It is apparent that these stars reveal similar sur- 
face enrichments, despite their differences in spectral types 
(Galactic: WC5-8; LMC: WC4) and stellar luminosities. 
The shift to (observed) higher luminosity is reasonably 
well reproduced in the non-rotating evolutionary models, 
although the shift in the (sole) rotating track is dramatic. 
Therefore, we must await results for a wider initial range 
of masses and rotational velocities before we can have con- 
fidence in the validity of these predictions. 
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Fig. 11. Comparison between the (C-|-0)/He values of 
WC stars in the LMC and Galaxy versus stellar lumi- 
nosity, together with evolutionary predictions for non- 
rotating stars (Meynet et al. 1994) for the two metallic- 
ities, plus results from an initially rotating 40Mo solar 
metallicity model from Meynet & Maeder (2000). Symbols 
(and BAT/WR numbers) are as for the previous figure. 



5.3. Differences in WC spectral types between Galactic 
and LMC stars? 

While we have established that WC stars in the LMC are 
typically more luminous, and possess lower mass-loss rates 
than Galactic counterparts, what is the origin of their 
systematically earlier spectral types,? This has been as- 
sumed to be due either to temperature differences and/or 
an abundance effect. 

Smith & Maeder (1991) adopted the hypothesis that 
the difference between the WC subtypes was primar- 
ily governed by a variable surface abundance ratio, 
(C-|-0)/IIe. This was based on recombination line studies 
of a sample of Galactic WC4-9 stars. However, only a sin- 
gle WC4 star was included in their sample, whilst recombi- 
nation studies of late-type WC stars are fraught with dif- 
ficulties (Hillier 1989). The present results reveal further 
that the expected (C-hO)/He ratios for WC4 stars (0.7- 
1.0) are not supported here (0.1-0.4), with all Galactic 
and LMC WC stars poorer in carbon and oxygen than 



12 



Crowther et al.; Study of LMC WC4 stars 



standard (non-rotating) evolutionary models would sug- 
gest. 

With regard to oxygen abundances, evolutionary pre- 
dictions for non-rotating massive stars indicate 0/0=0.1- 
0.22 by number at the stage in which 0.1< 0/He <0.3, 
in good agreement with 0/0=0.1-0.25 determined here. 
Predictions from evolutionary models in which rotation is 
considered (albeit for solar abundances only) indicate a 
fairly similar range (0/0=0.2-0.3 by number). 

Orowther (1999) proposed that the variety of early- 
type WO subtypes results principally from differing wind 
densities, rather than surface abundances or ionization. 
This seems to be supported in Fig. |l^. We have further 
illustrated this effect in Fig. ^ where we present UV and 
optical observations of HD 156385 (W07, Oalactic) and 
HD 32257 (W04, LMO). Olearly these distinct spectral 
types do not differ dramatically in their optical spectrum, 
other than the strength of the ill classification diagnos- 
tic at A5696. Hillier (1989) first emphasised the sensitivity 
of this line to stellar parameters, owing to the alternative 
decay mechanism from the upper level to 574A with a 
branching ratio of ~150. Oonsequently, A5696 will only be 
strong when A574 has a large optical depth. Note that 
other Oiii lines, for instance the feature at A6740, do 
not differ significantly between the two stars. Oonti et al. 
(1990) has previously remarked upon how other strong 
Oiii lines (A4650, A9710) do not closely follow the classi- 
fication Oiii line. In addition, the Hen A5412-Oiv A5471 
blend is remarkably similar for these two stars, such that 
their 0/He abundances are essentially identical. 

We include synthetic spectra for HD 32257 (from the 
present work) in Fig. |l^, plus two further models in which 
changes to the temperature and mass-loss rate are made. 
The temperature is reduced from 85kK to 80kK, and then 
to 75kK, whilst the mass- loss rate is increased from 1.5 
xlO^^Afo yr^i, to 2.1 xlO^^Af© yr"i, and then to 2.8 
xlO~^M0 yr~^, in order to closely match the parameters 
determined by Dessart et al. (2000) for HD 156385. From 
the figure, it is apparent that this small shift in parameters 
is sufficient for the spectral type to shift from W04 to 
W07, whilst other Oiii lines remain remarkably stable. 
Stellar temperature certainly plays a role - fixing the mass- 
loss rate and increasing the luminosity at constant radius 
also reduces the strength of A5696. In the UV, the effect 
of varying wind density is also evident. For example, the 
higher mass-loss rate of HD 156385 dramatically weakens 
OlV AA1548-51 through iron-group blanketing. 

From Fig. |l^ one might presume that mass-loss plays 
a principal role in determining spectral types for all WO 
stars, via the effect on Oiii A5696, However this is not al- 
ways the case. It has been established for several decades 
that early-type WO stars possess broader lines than late- 
type WO stars. Whatever its physical origin, this also has 
the effect of weakening the classification ill A5696 line. 
From recent, mostly unpublished calculations, it is ap- 
parent that W08 and especially W09 stars do have sig- 
nificantly lower stellar temperatures than other WO sub- 
types. Wind density does appear to play the principal role 



amongst W04-7 stars, because of the sensitivity of the 
classification ill A5696 line to mass-loss rate. Therefore, 
the spectral type amongst W04-7 stars, which we consider 
to represent 'early-type' WO stars is principally defined by 
wind density, with temperature a secondary effect and no 
role played by carbon abundance. W04 and WO stars def- 
initely possess weaker winds than other WOE stars, and 
probably have the highest stellar temperatures. Most likely, 
weaker winds and higher temperatures distinguishes WO 
stars from W04 subtypes. 

In comparison, the classification criteria for WN stars 
are much stronger indicators of ionization than those of 
WO stars. The sole exception is that Niv A4058/Niii 
A4634-41 is sensitive to nitrogen abundance at fixed tem- 
perature, as illustrated by Orowther (2000). This ratio 
now forms the basis for classification of early O stars 
(Walborn et al. 2002a). 

6. Wolf-Rayet Masses 

6.1. Present and final masses of WC stars 

Present LMO W04 masses can be calculated based on 
their stellar luminosities, according to Schaerer & Maeder 
(1992) or Langer (1989), revealing values in the range 12- 
18M0. However, these calculations neglect the fact that 
measured (photon) luminosities in WR stars will be less 
than that generated by nuclear processes in the deep in- 
terior, the so-called 'wind darkening' effect. Taking wind 
darkening into account results in a slightly larger mass 
(Heger & Langer 1996). These calculations are presented 
in Table 0, together with masses derived from Langer 
(1989). Note that the differences are much smaller than 
those reported previously due to the reduced mass-loss 
rates for WR stars once clumping is taken into account. 
We have confidence in our derived masses due to the close 
agreement between the mass of the W08 component in 
7 Vel obtained here from a previous spectral analysis 
{^9Mq) with that obtained from an orbital study of this 
binary system (9.5M0, De Marco et al. 2000). 

Further, in order to estimate final masses we make the 
following assumptions: 

— Lifetimes remaining prior to a supernova explosion can 
be adopted from recent evolutionary models. For ex- 
ample, HD 32125 (for which O/He=0.18 by number) is 
predicted to end its life in 3.7x10^ yr according to the 
O.AZq tracks from Meynet et al. (1994) irrespective of 
its initial mass (in the range 6O-I2OM0). Alternatively, 
Meynet & Maeder (2000) have presented evolutionary 
tracks for rotating (solar metallicity) massive stars - 
following their predictions at 4OM0 a much lower re- 
maining age is obtained - 0.9x10^ yr. Since the de- 
pendence of metallicity and initial mass remains to be 
established for rotating models we consider both pos- 
sibilities in Table 0. 

— We adopt a mass (luminosity) dependent mass-loss 
rate for the remainder of the star's life, scaled ac- 
cording to Eqn. 1 and using the Langer (1989) mass- 
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Table 7. Determinations of masses for Galactic and LMC 
WC stars using the relationship from Langer (1989) for 
stellar (photon) luminosities (MLsg)- Correction for wind 
blanketing effects discussed by Heger & Langer (1996) 
imply slightly higher (nuclear) luminosities (I/nuc) and 
masses (Mnuc)- We also provide predicted final prc-SN 
masses (Mpro-SN) based on estimates of lifetimes (r) re- 
maining, using alternative evolutionary models - Meynet 
ct al. (1994, OAZq and I.OZq) or Mcynct & Maeder (2000, 
I.OZq) as discussed in the text, plus the time dependent 
mass-loss rates scaled according to Eqn. 1. Stellar prop- 
erties for Galactic stars arc taken from Hillier & Miller 
(1999), Dessart et al. (2000), De Marco et al. (2000), 
Smartt et al. (2001). 
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luminosity relationship. For the above example with 
standard evolutionary tracks, the mass-loss rate is re- 
duced from 1.6x10"^ to 0.7x10"^ Mq yr~^ over the 
next 3.7x 10^ yr, ending its life with a mass of 11. IM©. 

Consequently, we consider final, prc-SN masses of 11 
19M0 for the entire LMC sample. Given the systemat- 
ically less luminous Galactic stars, we determine corre- 
spondingly lower final masses of 7 14M0 for these stars. 

In addition, we are able to estimate the total amount 
of carbon released by the star throughout its WC life- 
time which we have calculated as follows. For each epoch 
we have determined the mass fraction of carbon, Xq, via 



quadratic or linear fits to predictions of Xc versus remain- 
ing lifetime from Meynet et al. (1994) or Meynet & Maeder 
(2000). The amount of carbon liberated is then simply the 
product of the age interval, age dependent mass-loss rate 
(see above) and carbon mass fraction. In the above ex- 
ample, HD 32125 is predicted to lose a total of 2.2M0 of 
carbon during its WC phase using figures from Meynet et 
al. (1994) or 0.6Mq from Meynet & Maeder (2000). Our 
six WC4 targets are predicted to provide the LMC inter- 
stellar medium with 5-20 Mq of carbon over their WC 
lifetimes. 



6.2. What masses for Type Ic SN? 

It is possible to determine lower limits to the CO masses 
of Type Ic SN from modelling their light curves and fit- 
ting expanding models to their early spectra. For ex- 
ample, Iwamoto et al. (1998, 2000) present fits to three 
SN Ic's, revealing CO core masses in the range 2.IM0 
(SN 19941), IOM0 (SN1997ef), I3.8M0 (SN1998bw). 
Recently, Smartt et al. (2002) have discussed potential 
precursors to the type Ic SN 2002ap in M74, which in- 
clude a WC star. From an examination of pre-SN images 
of this galaxy (with an LMC like metallicity) , they de- 
duced an upper limit of Mb = — 6.3±0.5 mag for the im- 
mediate pregenitor. The brightest of our six WC stars has 
Mb = —5.3 mag, which corresponds to Mb = —5.8 mag in 
the Johnson system. None of these stars would have been 
visible on pre-SN images, adding substance to a possible 
WC precursor to SN 2002ap. 

The values presented above represent minimum values 
of the final progenitor mass, since an unknown fraction 
may have imploded to form a black hole or other compact 
object. Nevertheless, the CO ejecta mass of SN 1998bw fits 
in remarkably well with the predicted final masses of LMC 
WC stars determined here. Evolutionary predictions, plus 
model analyses, imply that He represents a non-negligible 
fraction of the final surface abundance of very massive 
stars. This is consistent with the initial He-free Type- 
Ic classification for SN 1998bw, since subsequent observa- 
tions did reveal the presence of Hei 1.083/im (Patat et al. 
2001). 

SN 1998bw was also remarkable in that it appeared to 

be CO- incident with GRB 980425. This Gamma Ray Burst 
was relatively sub-luminous, and it has been argued by 
Woosley et al. (1999) that such outbursts are atypical for 
GRBs. Nevertheless, SN 1998bw/GRB 980425 most likely 
corresponds to the collapsed core of a very massive star 
at the end of its WC phase. Recent XMM-Newton obser- 
vations of GRB 011211 by Reeves et al. (2002) provides 
further support for such an origin due to enriched Mg, Si, 
S, Ar, Ca and Ni in the X-ray afterglow. 

The lower CO mass of SN 1997cf is consistent with that 
of a luminous WC stars in our Galaxy, such as HD 165763 
(WRlll, WC5) for which we obtain a final mass of 10- 
12M0. The exceptionally low CO mass of SN 19941 has 
no analogue amongst Galactic or LMC massive WC stars. 
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unless most of the final CO mass imploded as a black 
hole/neutron star. 

7. Conclusions 

We have used far-UV, UV, optical and near-IR spec- 
troscopy together with line blanketed model atmospheres 
to re-evaluate the stellar parameters of six LMC WC4 
stars. We derive stellar luminosities which are a factor 
of ^2 higher than previously established for these stars 
(Grafener et al. 1998), and systematically higher than 
Galactic counterparts at known distance. Mass-loss rates 
are lower than previously determined due to the inclusion 
of clumping in our models. Several lines of evidence, in 
both WN and WC stars, strongly suggest that the winds 
of WR are strongly clumped. 

We find that the LMC WC4 stars possess very simi- 
lar surface abundances to Galactic WC5-8 stars, although 
their wind densities are systematically lower, by ^--^0.2 dex. 
We propose that the lower heavy metal content of the 
LMC is responsible for their lower mass-loss rates via 
their (still unproven) radiatively driven winds. A metal- 
licity dependence of ~ Z'^-^ would produce this weak ef- 
fect, comparable to that predicted for O stars as a func- 
tion of metaUicity {M oc Z°-^-° '^: Kudritzki & Puis 2000; 
Vink et al. 2001). This relatively minor difference quanti- 
tatively explains their different subtype distribution (re- 
call Fig.|l2|). Temperature appears to be a secondary cri- 
teria in distinguishing WC4~WC7 stars. Late-type WC 
stars, especially WC9, do appear to be rather different, 
with systematically lower stellar temperatures. Looking 
at the broader picture of the subtype distribution in other 
Local Group galaxies, other factors clearly come into play. 

If metaUicity were the sole factor which determined the 
wind strengths of WC stars, those galaxies with metallic- 
ities close to Galactic (e.g. M33 and M31) would have a 
large, late-type WC population, which is apparently not 
the case (Massey & Johnson 1998), unless these are all 
obscured in the visible by circumstellar dust. Recent ev- 
idence suggests that most (or all) WC9 stars are close 
binaries with time-dependent dust production from inter- 
actions between the winds of the WC and OB companion 
(e.g. Tuthill et al. 1999). Future quantitative studies of the 
properties of these and more distant WR stars will hope- 
fully help address these issues. Regardless of these details, 
the identification of a heavy metaUicity dependence of WR 
winds has a major impact on their relative contribution 
to the hard ionizing photons in young starburst galaxies 
(Smith et al. 2002). 

The association between spectral type and wind den- 
sity amongst LMC and Galactic WCE stars also has rele- 
vance to the long standing debate amongst [WC]-type cen- 
tral stars of Planetary Nebulae (CSPN). It has long been 
recognised that [WC] stars positively avoid WC5-7 spec- 
tral types, and either cluster around WC4/W0 or WC8- 
10 (e.g. Fig. 2 in Acker et al. 1996; Crowther et al. 1998). 
A wind density origin for the differences between WC4- 
7 stars (rather than temperature) can therefore explain 



the observed clustering as due to reduced wind densities 
amongst [WC]-type CSPN relative to massive Galactic 
WC stars, without the need to infer rapid evolution, for 
instance. Clearly, one still would need to understand why 
[WC]-type CSPN possess weaker winds than their massive 
cousins, but the differences would not need to be great - 
a factor of ^2 would suffice. 

Finally, we determine current masses of LMC and 
Galactic WC stars including the (minor) effect of wind 
darkening (Heger & Langer 1996). Based on remaining 
lifetimes from evolutionary predictions we estimate final 
pre-SN masses of ll-19Mo for LMC WC stars, and 7- 
14Mo for Galactic WC stars with known distances. These 
values are consistent with WC stars being the immediate 
precursors to luminous Type-Ic supernova explosions, in- 
cluding SN 1998bw and SN 1997ef. 
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Fig. 4. Comparison between the de-reddened (i;B-v=0.07(Galactic)-|-0.04(LMC)) flux distribution of HD 37026 (solid, 
black) and theoretical predictions (dotted, red) for = 90kK, log L/Lq=5.65, logM/(MQ yr~i) = -4.5, i;oo=2900 
kms~^, C/He=0.32, O/He=0.05 by number. For the FUSE region, we use an atomic hydrogen column density of 
log Ai"(H/cm~^) = 21.0, whilst correction for molecular hydrogen is taken from H2 models by Tumlinson et al. (2002) 
using a Doppler width (6 = 5.4 kms~^) and total column of log (H2/cm~^) = 15.45, considering levels J=0 to 3. 
Ordinate units are ergcm"^ s~^ , whilst abscissa units are A except for jira in the bottom panel. 
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Fig. 6. Comparison between de-reddened observations of HD 37026 (solid, black) and synthetic spectra for a variety 
of oxygen abundances - 5% by mass (dotted, green), 10% (dashed, red) and 20% (dot-dash, blue). OlV AA3403-13 
is rather insensitive to abundance, whilst the strong A2900 feature is a blend of CiV A2906, 2918 and Oiv A2921. 
Ordinate units are ergcm^^ s^^ A~^. 
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Fig. 12. Rectified optical (upper panel) and UV (lower panel) spectroscopy of HD156385=WR90 (WC7, Galactic) and 
HD32257=Br8 (WC4, LMC), together with synthetic spectra from a series of three models, with log {L/Lq)=5A9, 
Uoo=2300 kms~^, C/He=0.35 and O/He=0.06, except that (bottom to top) the temperature:mass-loss rate is 
75kK:2.8xlO-5 Moyr"!, 80kK:2.1xlO-5 Mgyr"!, and 85kK:1.5x IQ-^Mgyr-i, respectively. Remarkably, these small 
changes in stellar parameters closely match the optical and UV spectra of these WC4 and WC7 stars, with C iii A5696 
most sensitive to small changes in physical parameter (compare A5696 with A6740 in the upper panel). 



